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ABSTRACT 

The prompt emission of gamma-ray bursts probably comes from a highly relativistic 
wind which converts part of its kinetic energy into radiation via the formation of shocks 
within the wind itself. Such " internal shocks" can occur if the wind is generated with 
a highly non uniform distribution of the Lorentz factor. We estimate the expected 
photosphcric emission of such a relativistic wind when it becomes transparent. We 
compare this thermal emission (temporal profile + spectrum) to the non-thermal 
emission produced by the internal shocks. In most cases, we predict a rather bright 
thermal emission that should already have been detected. This favors acceleration 
mechanisms for the wind where the initial energy input is under magnetic rather than 
thermal form. Such scenarios can produce thermal X-ray precursors comparable to 
those observed by GINGA and WATCH/GRANAT. 

Key words: Gamma-rays: bursts - Radiation mechanisms: non-thermal - Radiation 
mechanisms: thermal - Hydrodynamics - Relativity 



1 INTRODUCTION 

The cosmological origin of long duration gamma-ray 
bursts (hereafter GRBs) has been firmly established s ince 
the discovery of th eir optical counterparts in 1997 (van 
Paradijs et al. 1997). These late and fading counterparts, 
the so called afterglows, have now been detected in many 
bursts, and in different spectral ranges : X-rays, optical 
and radio bands. The redshift has been mesured for about 



very similar systems made of a stellar mass black hole sur- 
rounded by a thick torus. The collapsar model seems to be 
favored in the case of long bursts by observational evidences 
that GRBs are located well inside their host galaxy and 
often associated to st ar-forming regions ( Paczynski 1995 ; 
Djorgovski et al. 2001). The released energy is first injected 



into an optically thick wind, which is accelerated via an 
unknown mechan is m, probably involving MHD processes 



20 GR Bs from z 



0.43 Lo z 



4.5. The corresponding 



isotropi c equivalent energy radiated by these GRBs in the 

gamma-ray range goes from 5 10 51 erg to 2 10 54 erg. The 



(Thompson 1994; Meszaros & Rees 1997; Spruit et al 



2001) and becomes eventually relativistic. The existence 
of such a relativistic wind has been directly inferred from 



beaming facLoi thai has Lo be taken into account Lo obtain 
the rea l amount of radiated energy can be deduced from 

afterglow obse rvations (achromatic break in the lightcurve, 
Rhoads (1999)). Current estimates lead to a total energy 
radiated in gamma-rays of about 0.5 — 1 10 51 erg (Frail et 
al. 2001). The most discussed scenario to explain the GRB 
phenomenon is made of three steps : 

Central engine : The source of GRBs must be able to 
release a very large amount of energy in a few seconds. 
The two most popular candidates are either the merger of 
compact objects (neu tron star binaries or neutron star 



the observa tions of radio scintillation in GRB 970508 (Frail 



et al. 1997) and is also needed to solve the compactness 



problem and avoid photon-photon annihilation along the 
line of si ght. Average Lorentz fac t ors larger than 100 ar e 
required (|Baring fc Harding 1997|; [Lithwick fe Sari 200l[). 



The next two steps explain how the kinetic energy of 
this relativistic wind is converted into radiation at large 
distances from the source, when the wind has become 
optically thin. 

Internal shocks : the production of gamma-rays is usually 
associated to the formation of shocks within the wind itself 



black hole systems ( Narayan et al. 1992 



Mochkovitch 



et al. 1993)) or the gravitational collapse of a massive star 



into a black ho le (collapsars/hypernovae ( Woosley 1993 



Paczynski 1998)). Such events lead to the formation of 



(Rees & Meszaros 1994). Such internal shocks can appear if 
the initial distribution of the Lorentz factor is highly vari- 
able, which is very likel y considering the unsteady na ture of 
the envisaged sources (MacFadyen fc Woosley 1999). This 
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model has been studied in d e tails (Kobayashi et al. 1997 
Daigne fc Mochkovitch 1998|, |2000D. The main difficulties 



which are encountered are a rather low efficiency for the 
conversion of the wind kinetic energy into gamma-rays 
(a few percents only) and problems in reproducing with 
synchrotron e mission the slope of the low energy part of 
the spectrum (GJhisellini et al. 2000). Despite this difficulty, 



the model can successfully reproduce the main features of 
the bursts observed by BATSE. 

External shock : the relativistic wind is decelerated 
later by the external medium. This phase of deceleration 
is probably the best understood of the thr ee steps and 



repr odu ces very well the afterglow properties (Wijersetal 
1997|)~lrhe 

dynamics of the wind during the deceleration 
phase is d escribed by the solution of t he relativistic Sedov 
problem ( |Blandford fc McKee 197t ) and the observed 



afterglow is due to synchrotron emission produced by 
relativistic electrons accelerated behind th e strong forwar d 
shock propagating in the external medium (Sari et al. 1998). 



The work presented in this paper focuses on the prompt 
emission. The spectrum of this emission as observed by 
BATSE and Beppo-SAX is non-thermal and is w ell fitted 
by the 4-parameter "GRB-function" proposed by Band et 
al. (199^). This function is made of two smoothly connected 
power-laws. This non-thermal emission probably originates 
from the radiation of a population of highly relativistic 
electrons accelerated behind the shock waves propagating 
within the wind during the internal shock phase. 



r.«(t ln ,;t) 



r(t, „,';t+At) 



First layer 



Layer emitted 
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Wind production 

Figure 1. Sketch of the plane t-r : the solid line arrows show the 
paths of the layers produced by the source at radius r acc from 
t = to t = tyf. The dashed line arrow shows the path of the 
photons emitted at time t by the layer produced at time f; n j. 
These photons will cross the layer produced at time . < t ln j 
after a duration At. They escape when they cross the first layer 
of the wind produced at t = 0. 



Prior to the internal shock phase, the relativistic wind 
has to become transparent. At this transition, a thermal 
emission is produced, that could contribute to the observed 
prompt emission. Parts of the wind can also become 
opaque at larger radii if internal shocks create pairs in 
large number. These opaque regions can produce additional 
t hermal components wh en they become transparent again 
(Meszaros & Rees 2000). Other thermal contributions can 
be expected, for example when the jet breaks out at the 
boundary of the stellar en velope in the collapsar scenario 
( Ramirez- Ruiz et al. 2002). In this paper, we restrict our 
analysis to the photospheric t hermal component. A sim ilar 
problem has been studied by Lyutikov & Usov (200C) in 
the different context of strongly magnetized winds emitted 
by rapidly rotating pulsars. 

The paper is organized as follows : in sect. |^ we obtain 
the position of the photosphere of a relativistic wind with a 
highly variable initial distribution of the Lorentz factor, as 
expected in the internal shock model. We then compute the 
corresponding photospheric thermal emission in sect. ^ and 
compare it to the non-thermal emission from the internal 
shocks in sect. The results are discussed in sect. ^ and the 
conclusions are summarized in sect. H. 



2 THE PHOTOSPHERE OF A RELATIVISTIC 
WIND 

2.1 Photospheric radius 

We do not discuss in this paper the nature of the source 
which is initially responsible for the energy release leading 



to the gamma-ray burst. We suppose that a relativistic 
wind carrying the energy has emerged from the source, 
with an average Lorentz factor T > 100. We assume that 
the acceleration is complete at a distance r acc from the 
source where the ultra-relativistic wind is characterized by 
an energy injection rate _E(tj n j) and an initial distribution 
of Lorentz factor r(ti n j). This corresponds to a mass flux 
M(tinj) = E(t in j)/r(t in j)c 2 , with E and M being the 
isotropic equivalent energy and mass injection rates. This 
wind production process lasts from t[ n j — to tinj — t w (all 
these quantities are defined in the fixed frame of the source) . 

In this section we are interested in computing when the 
layer emitted by the source at ti n j will become transparent. 
We assume that the wind is still optically thick at r acc , and 
that it becomes transparent before the internal shock phase 
and before it is decelerated by the external medium. We can 
then consider that each layer is evolving with a constant 
Lorentz factor so that at time t, the layer emitted at ti n j is 
located at 



r(tinj,t) 



C(t - t inj ) 



(1) 



Let us consider photons emitted at t by the layer ejected 
by the source at tinj (see fig. |l|). If they escape from the 
relativistic wind, these photons will have to cross all the 
layers emitted from t[ n j = to t[ a ; = ti n j. Precisely, they 
cross the layer ejected at t[ n j after 



Ar 



At~2r 2 (t; nj )— , 



(2) 



where Ar is the spatial separation at the emission time t 
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between the layer produced at t[ n j and the emitting layer 
produced at time tinj . This distance is the initial separation 
c (tinj — ti n j) P ms a correction growing with time due to the 
difference of Lorentz factor between the two layers : 



Ar 



" (t' ioi ,t) - r (t in j,t) 




i r 2 (t; nj )-r 2 (t mj ) 
2 r=»(c)ra(fcj) c( 



tin 



(3) 



The first term is very close to the initial separation and 
the second term is small as long as the process we consider 
takes place well before the internal shock phase. The photons 
escape from the wind when they cross the first layer emitted 



at time 



at t' iB 



t esc ~t + 21^(0) t m + 



i r 2 (o)-r 2 (t inj ) 

2 r2(0)T2 (i iaj ) 



(t-tim)) (4) 



and at radius r csc (ti n j;t) — r(ti n j;t) + c(t csc — t). The cor- 
responding distance is 2r 2 (0) times larger than the initial 
separation between the emitting layer and the front of the 
wind. The total optical depth for these photons is given by 

rr osc (t in j,t) 

r(t inj ,t) = / dr(r) . (5) 

The elementary cont ribution dr(r) to the op tical depth is 
a Lorentz invariant ( Abramowicz et al. 1991 ) and is more 
easily estimated in the comoving frame of the layer crossed 
by the photons at r : 



dr(r) = Kp dZ'(r) 



(6) 



where « and p' are the opacity and the comoving density 
of the layer. The length dl'(r) is computed by a Lorentz 
transformation from the fixed frame to the comoving frame 
of the layer. We take into account the fact that when photons 
cover a distance dr, the corresponding duration is dt = dr/c 
so that 

1 M , dr 
~ 2T 2 

The comoving density is given by 



dZ'(r) = r(dr-vdt) 



dr : 



2r 



(7) 



1 M 



P - T. 



r 47rr 2 c 



(8) 



Here all the physical quantities like V, M, etc. have the value 
corresponding to the layer crossed by the photons at r, i.e. 
the layer emitted at t( n j solution of (from eq. ^ 



r = r(t inj ,t) + 2r 2 l y t' ini )Ar , 



(9) 



where Ar is given by eq. |^. The final expression for the total 
optical depth is 



r(t in j,t) 



: (tin.pt) 



KM 



87rr 2 r 2 c 



dr 



(10) 



We define the photospheric radius r pn (tinj) of the layer 
emitted at ijnj by 

r p h (tinj) = r(ti n j,f) with r(t in j,t) = 1 . (11) 
To estimate this radius, we still need to specify the opacity. 
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Figure 2. A single pulse burst : initial di strib ution of the 
Lorentz factor (example considered in sec. 2.4). The Lorentz 



factor r is plotted as a function of the injection time fi n j. As a 
constant energy injection rate E is assumed, the masses of the 
fast (r = 400) and "slow" parts (r = 100 -> 400) are equal. 



We consider here the phase when the acceleration is com- 
plete. The internal energy has already been almost entirely 
converted into kinetic energy. Pairs have annihilated and do 
not contribute to the opacity. Then the optical depth is due 
to the ambient electrons and the opacity is given by the 
Thomson opacity k = kt. In the following, when a numeri- 
cal value is needed, we use k = 0.2 (i.e a number of electrons 
per nucleon Y e — 0.5). 



2.2 The case of a constant Lorentz factor 

In the case of an homogeneous wind where E, M and T are 
constant, we have 



KMtin 



ntinj,t) ~ 47rr ( tin . )t ) resc ( tiaj)t ) 

There are two limiting cases : 



T(tinj,t) 



KMtin 



Anr 2 (tinj,t) 



if r(tinj,t) > 2r 2 ctin 



(12) 



(13) 



T ^ * SnJr^t) if r ^ « 2r2 ^ ■ (14) 

The corresponding photospheric radius of the layer ejected 
at tinj is given by 



, KMtinj A' 7 ' £ 167TC 4 

r ph (tini) I — ; I if ^r: — > 



4tt 



r 5 t;n 



r p h (tin 



KM . r E 16ttc 4 

if ^T— < 



87rcr 2 r 5 t ; „ 



(15) 



(16) 
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where we have replaced M by E/Fc 2 . The condition 
E/F 5 t inj < 16ttc 4 /k reads 



E 52 200 

F^inj ^0.2 



(17) 



which is usually true. Here -E52, and K0.2 are respectively 
E, F and « in unit of 10 52 erg/s, 10 2 and 0.2. Then, the 
photospheric radius is the same for all the layers and is given 
by (eq. @) : 

(18) 



r ph =3.0 10 12 k,o. 2^52^2 3 cm 



If we estimate r acc by the saturation radius Fro which is 
predicted in the fireball model, we get 



9 10 s nir 2 cm 



(19) 



for a typical initial radius ro taken to be the last stable orbit 
at three Schwarzschild radii around a non rotating black hole 
of mass Mbh = 10/iiMQ. It is clear that r p h is much larger 
than r acc as expected. 

2.3 The case of a variable Lorentz factor 

We now consider the case where the initial distribution of the 
Lorentz factor is v ariable. We use the simple mod el that has 
been developed by Daigne & Mochkovitch (1998). The wind 
is made of a collection of "solid" layers ejected regularly on a 
time scale At[ n j with a Lorentz factor, a mass and an energy 
Fi, Mi — MiAti n j and Ei — EiAtinj where i = 1 corresponds 
to the first layer produced at ti n j = 0. Photons emitted by 
the layer io when it is located at r j travel through a total 
optical depth 



, x _ K \ ~> Ei I J_ 1 

{ri °> ~ 8nc 3 2^ rf I rj n ~ r° ut 

i<in 



(20) 



been considered in Daigne & Mochkovitch (1998, 2000) and 
leads to a typical single pulse burst. We use Ati n j = 2 10 -3 s 
so that the wind is made of 5000 layers. For each layer i we 
compute the photospheric radius r p h i and the radius r csc ( 
where the photons emitted at the photosphere escape from 
the relativistic wind. The result is plotted in fig |§] as a func- 
tion of the mass coordinate Mi. Notice that except for the 
front of the wind, the approximate value of r p h given by 
eq. ^ works extremely well. An interesting result is that 
the deepest layers in the wind become transparent before 
the layers located at the front. This is due to the fact that 
photons emitted by these layers cross the front at larger 
radii when the density has already strongly decreased (one 
can see that r osc for these layers is larger than at the front). 
The photospheric radius goes from ~ 4.7 10 10 cm to ~ 
3.0 10 12 cm. We can check now that this is well before the 
internal shocks form or the deceleration of the wind by the 
external medium becomes efficient. The typical radius of the 
internal shocks is given by 



(24) 



where t var is the characteristic time scale for the variations of 
the Lorentz factor and / is a numerical factor depending on 
the details of the initial distribution of the Lorentz factor (/ 
will be smaller for high contrasts of F) . For a typical average 
Lorentz factor F > 100 we have 



ns 



3 10 14 / Fjtvar cm 



(25) 



and we immediately see that except for very small values 
of / or very short time scales i var , the typical radius of the 
internal shocks is larger than the photospheric radius. The 
deceleration of the wind by the external medium occurs even 
further away, except in very dense wind environments. 



where rf 1 and r° ut are the radii at which the photons enter 
and escape the layer i. We have (from eq. [)]) 



2Ff 



C (i - i) Atinj + 



2T1 



and 



r° ut =i C + 2F 2 C At in 



(21) 



(22) 



In the following, we use the exact formula ( |20[ ) to compute 
the optical depth and we solve numerically i"(r p h iQ ) = 1 to 
get the photospheric radius of the layer io. An approximate 
value is obtained under the assumption that the opacity is 
dominated by the contribution of the layer io where photons 
are emitted. We then have : 



K P 

approx nj±j Z q 



87rc a rl 



(23) 



for r ph iQ < 2F 2 ct inj (with t inj = (i - l)At inj ). This is 
usually true except for the first layers (tinj — * 0). 



2.4 Example : a single pulse burst 

We consider the case of a relativistic wind ejected dur- 
ing t w = 10 s with a constant energy injection rate E — 
10 52 erg/s and an initial distribution of Lorentz factor repre- 
sented in fig. 0. Such a simple initial distribution has already 



3 TIME PROFILE AND SPECTRUM OF THE 
PHOTOSPHERIC EMISSION 

3.1 Photospheric luminosity 

In the framework of the fireball model, the temperature and 
luminosity of a layer at its photo spheric radiu s r pn are given 
(in the fixed frame) by (see e.g. Piran (1999)) 

/ r . , \ -2/3 

kT ph ~ kT° ( ) (26) 
and 

L ph ^E(^Y 2/3 (27) 



(in this section we omit the index i. Everything applies to 
each layer). The radius r acc is the saturation radius defined 
by eq. |l9| and the initial b lackbody temperature of the layer 
is ( |Mesziros fc Rees 200o| ) 
-1/2 



^52 



MeV. 



(28) 



Deviations from the predictions of the standard fireball 
model are however possible. The central engine of gamma- 
ray bursts is still poorly understood and the acceleration 
mechanism not clearly identified. A large fraction of the 
energy released by the source may be for instanc e initially 
stored under magnetic form ( Spruit et al. 2001 ). In this 
case, the wind is not as hot as in the standard fireball 
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Figure 3. A single pulse burst: photosphere of the relativistic wind (example considered in sees. 2.4 and 3.3). Left: the 



photospheric radius r p h is plotted as a function of the mass coordinate of layers within the relativistic wind (solid line). The approximate 
value eq. ^ (dotted line) is also shown : the only noticeable difference is at the front of the wind (M — * 0). The dashed line shows the 
radius at which the photons emitted at the photospheric radius of each layer escape from the wind. Right : the photospheric luminosity 
L p h and observed temperature fcT p h /(l + z) are plotted as a function of the arrival time t a for A = 1. A redshift z = 1 is assumed. 



model and the photospheric luminosity is also smaller. An 
extreme case would be the magnetic acceleration of a cold 
wind where the photospheric temperature and luminosity 
are negligible. 

Whatever the physics of this early phase may be, it should 
necessarily have the two following properties in common 
with the standard fireball model : 

(i) The acceleration mechanism must have a good effi- 
ciency. The observed isotropic equivalent gamma-ray lumi- 
nosity L 1 is indeed very high. To account for it, the internal 
shock model requires a isotropic equivalent kinetic energy 
flux 

J"l 0.1 

where L 7 51 is the observed isotropic equivalent gamma-ray 
luminosity in unit of 10 51 erg/s and / 7 „ t is the efficiency for 
the conversion of kinetic energy into gamma-rays in unit of 
0.1. As this kinetic energy flux is already very high, we can- 
not expect the source to release much more energy. There- 
fore the need for an efficient acceleration is unavoidable. This 
means that beyond r acc , the energy flux is completely domi- 
nated by the kinetic energy flux, like beyond the saturation 
radius in the standard fireball model. The main difference 
may probably be the value of r acc compared to the standard 
saturation radius ~ IYo. 

(ii) Beyond r acc the wind experiences a phase of adiabatic 
cooling due to spherical expansion. An efficient acceleration 
indeed implies that the wind is still optically thick at r acc . 
In this case fcT p h and L p h decrease as r~ 2 ^ 3 beyond r acc 
like in the sandard fireball model. Only the initial value of 
the temperature and the internal energy density at the end 



of the acceleration phase can be different from those of the 
standard fireball model. 

To account for our poor knowledge of the physical process 
responsible for the acceleration of the wind, we define A as 
the fraction of the energy which is initially injected under 
internal energy form. In the standard fireball model A = 1 
whereas < A < 1 in the other possible cases. With this 
definition we have 

kT ph ~A 1/4 fcT°(^)~ 2/3 (30) 

V face / 

and 

L ph ^Xe(^l)~ 2/3 . (31) 

The acceleration radius r acc may differ from the saturation 
radius given by eq. [Hj. However, we will show below that 
the relevant quantity to estimate the photospheric emission 
is the ratio L p h / kT p h which does not depend on r&cc , as 
long as the photospheric radius is large compared to r acc . 

3.2 Spectrum, count rate and arrival time of the 
photospheric emission 

We suppose that the photosphere radiates as a blackbody at 
temperature fcT p h . This is clearly a simplifying assumption 
as scattering processes may play an important role when 
the opacity is t>1. However, we believe that the possible 
deviations from a pure blackbody will not change our main 
conclusions. We also neglect corrections in the spectrum due 
to angular effects affecting photons originating from differ- 
ent regions of the emitting shell. We then consider that the 
emitted photons have a Planck distribution which is in the 
source frame : 
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Figure 4. A single pulse burst : spectrum and time profile (example considered in sees. 2.4 and [j,3| ), A redshift z = 1 is assumed. 
Left : the spectrum E 2 n(E) is plotted as a function of the photon energy in keV. Right : the count rate is plotted as a function of arrival 
time in the 3.5—8.5 keV band (top) and in the 50—300 keV band (bottom). The dashed line corresponds to the photospheric emission, 
the dotted line to the non-thermal emission from the internal shocks and the solid line to the total emission. The global spectrum is 
dominated by the luminous contribution from the photosphere of the rapid part of the wind, peaking at about 3.92 kT p ^ ~ 370 keV. 



dn ph {E) 



E 



AEAt J Planck (fcTph ) 4 cxp (^L-)-l 



(32) 



where Tpi anc k = J Q e x P x-i dx = ^ . Taking into account 
the redshift z of the source, the observer will detect a photon 



flux at energy E (observer frame) which is given by 

4 



C ph (E) 



^Planck 47rDg V kT ph 



£ph [ 1 +z 



E- 



exp 



,(33) 



where Dl is the luminosity distance at redshift z. The cor- 
responding count rate in the energy band [Ei; E2] is 



°12 — 



L p h 1 + z T 1 



p ll 



ph -^-Planck 



(34) 



where If? = J* 2 cxp l_ 1 dx and Xi )2 = (1 + z)E ia / kT ph . 
It is interesting to notice that the ratio L p h / kT p h depends 
neither on the shell radius nor on the saturation radius : 



P h 



5.0 10 57 X 3/4 E% A £' 2 ph/s 



(35) 



Then the count rate in the energy band \E\ \ E2] is given by 
Cg ^ 4.0 i±i A 3 / 4 i? 5 3 / 4 ^ Ph/cm 2 /s , (36) 

-^-Planck 

where D28 is the luminosity distance Dl in unit of 10 28 cm. 
The emitted photons will be detected at the arrival time i a 
(relatively to a signal travelling at the speed of light) : 

r p h 

C 



ta = t 



ph 



(37) 



where £ p h is the time when the layer reaches the radius r p h . 
We get 



ta — ^ini ~t~ 



(38) 



mj ' 2T 2 c 

With the approximate value of r p h given by eq. |23L we have 



iinj + 



kE 
8tvc 4 T 5 



(39) 



We already checked in the previous section that the sec- 
ond term is negligible compared to tiny Then t a — tiny The 
spreading of arrival times over a duration At a — r p -^ /2cF 2 
due to the curvature of the emitting surface is also negligible 
for the same reason (of course these estimations of t & have 
to be multiplied by 1 + z in the observer frame to account 
for the redshift). The fact that i a — (1 + z)tinj shows that 
the time profile of the photospheric emission, if observed, 
would provide a detailed direct information about the ini- 
tial distribution of Lorentz factor in the wind. 



3.3 Example : a single pulse burst 

We consider the same distribution of the Lorentz factor 



and injected power as in sec. 2.4 and we now compute the 



thermal emission of the photosphere for a standard fireball 
(A = 1). Fig. ^| shows the luminosity and the temperature at 
the photosphere as a function of the arrival time of photons. 
We did not use the approximations given by Eqs. p6|-p7| 
which are strictly valid only for r p h 3> r aC c but we used 



the ex act solution of the fireball equations (see e.g. Piran 
(1999])). We adopted a redshift z — 1. Fig. shows the cor- 



responding integrated spectrum of the photospheric emis- 
sion and the time profile in two energy bands : 3.5-8.5 keV 
which is one of the X-ray bands of Beppo-SAX and 50- 
300 keV which is the 2+3 gamma-ray band of BATSE. The 
photospheric emission of the "slow" part (ti n j = — > 4 s and 
t a — (1 + z)t[ n j — — > 8 s) has a temperature kT p h increas- 
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Figure 5. The ratio of the count rate due to the photospheric emission over the count rate due to the internal shocks is plotted 
as a function of the peak energy of the non-thermal spectrum (observer frame) for different values of the photospheric temperature 
fcTph /(1 + z) in the observer frame. The spectrum of the non-thermal emission from the internal shocks is computed using the GRB- 
function with a = —1.0 and j3 = —2.25. The following parameters have been adopted : A = 1, E52 = 1, fj = 0.1, fii = 1. Left: X-ray 
band 3.5—8.5 keV; Right: gamma-ray band 50—300 keV. 



ing from 4.2 to 94 keV. It initially produces a pulse only- 
visible in the X-ray band (t a ~ — > 3 s). Then, the count 
rate rises in the gamma-ray band, reaches a maximum at 
~ 5 s when kT pa — 48 keV and starts to decrease (although 
the temperature is still increasing) because the peak energy 
(~ 3.92 fcTph ) becomes larger than 300 keV. The rapid 



part (ti n 



10 s and t a 



(1 + z)t ini = 8 -» 20 s) has 



a constant temperature of 94 keV, so that the count rate is 
constant and mainly visible in the gamma-ray range. 



4 COMPARISON WITH THE EMISSION 
FROM THE INTERNAL SHOCKS 

4.1 Time profile and spectrum of the emission 
from the internal shocks 

We now estimate the count rate due to the emission of 
the internal shocks. Their luminosity can be written as 
Lis — ft E, with / 7 ~ fdUe- The efficiency /d of the dis- 
sipation process is the fraction of kinetic energy which is 
converted into internal energy behind the shocks; a c is the 
fraction of the internal energy which is injected into relativis- 
tic electrons, which then radiate to produce the gamma-ray 
burst, with a radiative efficiency / ra d which is assumed to 
be very close to 1. We do not discuss here the details of the 
radiative processes and we simply assume that the emitted 
photons h ave a spectral dis tribution given by the "GRB- 
function" (|Band et al. 1993|) : 



dn IS (£Q 
dEdt 



B 



(40) 



The peak energy E p is defined as the maximum of 
E 2 dn is (E)/dE/dt and is measured here in the source frame, 
so that the peak energy in the observer frame is E p / (1 + z). 



The function B(x) has two parameters, the low and high 
energy slopes a and (5 and is given by 



B(x) = 



x a exp(-(2 + a)a;) if x < §=§ 

x<i (§T^) a ^ ex P(fl — a) otherwise 



(41) 



The integral Teand = J + °° xB(x)dx depends only on a and 
P which we assume to be constant during the whole burst. 
The observed photon flux at energy E is given by 

and the corresponding count rate in the energy band \E\ ; E2] 
is 



1 T z { E \ — ^ 

Clt = "Dl7 fl 01 «keV j T^t 
where lf 2 = B(x)dx with xi, 2 = (1 + z)E 1)2 /E p . 



ph/cm 2 /s,(43) 



4.2 Comparison with the photospheric emission 

We now define R12 as the ratio of the count rate due to the 
photospheric emission over the count rate due to the internal 
shocks : 



R _ 1 « \3/4 f -1 p-l/4 1/2 
H12 — l.o A / 7 Q1 E 52 n x 



E n 



T 7-ph 
-^Band -^12 

200 keV TpianekZj! 



(44) 



Fig. q shows the value of R12 for two energy bands (X- and 
gamma-rays) as a function of E p assuming different values of 
the photospheric temperature kT pa . We have adopted a = 
— 1.0 a nd — —2.25 which are the typical slopes observed in 
GRBs (Preece et al. 2000). It is clear that the photospheric 
emission will show up in a given band (i?i2>0.1) when 
the observed temperature kT pn /(1 + z) of the photosphere 
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Figure 6. A single pulse burs t wi th a less luminous photosphere : spectrum and time profile (example considered in sec. 1.3). 
The wind is the same as in sec. 2.4 but we now assume A = 0.01. The redshift is z = 1. Left : the spectrum E 2 n(E) is plotted as a 
function of the photon energy in keV. Right : the count rate is plotted as a function of arrival time in the 3.5-8.5 keV band (top) and 
in the 50-300 keV band (bottom). As in fig. |l| the dashed line is the contribution of the photospheric emission, the dotted line is the 
contribution of the internal shocks and the solid line is the total count rate (photosphere + internal shocks). 



crosses this band. These results show that with the pre- 
diction of the standard fireball model for the photospheric 
temperature and luminosity, it is very difficult to prevent the 
photospheric emission from being easily detectable^ either 
in the X-ray or gamma-ray range. The presence of a bright 
thermal component is not supported by the observations : 
the gamma-ray burst prompt emission, as seen by BATSE, 
is clearly non-thermal. Concerning the X-ray emission, espe- 
cially at the beginning of the burst, more observations with 
better spectroscopic capabilities than Beppo-SAX will be 
necessary to check wether a thermal component is present 
or not. 



range, eq. |44| indicates that either / 7 must increase or 
A must decrease. The first solution is then to have more 
efficient internal shocks. In the example presented in fig. Q 
the efficiency is low : / 7 ~ 0.02. As there are many 
uncertainities in the radiative processes leading to the 
observed gamma-ray emission, one could hope that f~, 
is indeed very close to /d, the fraction of kinetic energy 
dissipated in internal shocks. However, the efficiency /d will 
never exceed a few 10~ . Therefore, even in the ideal case 
where / 7 ~ /d, one cannot expect to have / 7 larger than 
0.3-0.4. We have checked that this is not enough to avoid a 
easily detectable photospheric emission. 



4.3 Example : a single pulse burst 

We have computed the emission of the i nte rnal shoc ks in 
the single pulse burst considered i n sees. 2.4 and 3.3 using 



the si mple model developped by Daigne & Mochkovitch 
(1998). The result is shown in fig. W. The equipartition 



parameters have been chosen so that the peak energy of 
the emission from the internal shocks is 200 keV. As the 
photosphere reaches a temperature kT p ^ /(l + z) as high as 
94 keV, the study made in the previous subsection predicts 
that the photospheric emission should be easily detectable, 
which is clearly visible in fig. ^. 

To recover a burst which is dominated by the non- 
thermal emission of the internal shocks in the gamma-ray 



The only remaining solution is then to assume that the 
photosphere is less hot and luminous than what is predicted 
in the standard fireball model, i.e. to decrease the value of 
A in eqs. |3o|-|3l|. To have a spectrum entirely dominated 
by the non-thermal component we have to adopt A < 0.01, 
which means that less than 1 percent of the energy initially 
released by the source is injected under internal energy 
form into a "standard" fireball. Such a situation is shown in 
fig. |^ where we have computed the photospheric emission 
of the single pulse burst with A = 0.01. It is impossible 
to detect the thermal emission of the photosphere neither 
in the gamma-ray profile nor in the global spectrum. 
However, in the X-ray range, one can notice that during 2 
seconds before the internal shock emission starts, there is 
a thermal precursor whose intensity is about 8 percents of 
the intensity at maximum in the main pulse. 



Notice that E p is defined as the peak energy of the non-thermal 
emission of the internal shocks and is of course no more the peak 
energy of the total observed spectrum when the photospheric 
emission is dominant. 



We have finally considered the effect of pair creation 
during the internal s hock phase. The optica l depth for pair 
creation is given by ( Meszaros fc Rees 200C ) : 
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U± Lis 0~T 

4nr(m e c 2 )T 3 c 



(45) 



where a± is the fraction of the energy radiated in photons 
above the pair creation threshold. Our internal shock model 
allows the computation of t± at each shock radius. For the 
example considered here, t± never exceeds 6 10 -2 . For larger 
E and / or smaller Lorentz factors (in this case, if E is in- 
creased by a factor of ~ 20 or if all Lorentz factors are 
divided by ~ 2), r± increases and pair creation can become 
important, especially for shocks occuring at small radii. We 
do not compute in this paper the detailed internal shock 
spectrum for this case where an additional thermal compo- 
nent can be expected. 

5 DISCUSSION 

5.1 X-ray thermal precursors 

In the GRBs observed by BeppoSAX the X-ray and 
gamma-ray emission usually start simultaneously or the 
gamma-ray emission starts earlier. Usually, no ev idence is 
found for a thermal component in the spectrum (Frontera 
et al. 2000). Then, the prompt X-ray emission is probably 
due to the internal shocks like in the gamma-ray range. This 
implies that the photospheric emission must be present in 
these bursts only at a very low level, i.e. A 1 as explained 
in the previous section. However, in at least one case - 
GRB 990712 - evidence was found in the spectrum for the 
presence during the b urst of a weak therm al component of 
temperature 1.3 keV (Frontera et al. 2001). In complement, 



a X-ray pre cursor activity has bee n detected in a few GRBs 
by GINGA (iMurakami et al. 1991ft and WATCH/GRAN AT 



(Sazonov et al. 199^ In the observations carried out 



with the GRB detector onboard the GINGA satellite, 
X-ray precursors were detected between 1.5 and 10 keV 
in about one third of the GRBs. The spectrum of these 
X-ray precursors could be approximated by a black-body 
with temperatures between 1 and 2 keV. The WATCH 
catalog also includes several GRBs with X-ray precursors 
detected between 8 and 20 keV. As can be seen in the time 



profiles of these bursts (sazonov et al. 1998), the X-ray 
precursor usually has a duration which is about 20-50% 
of the duration of the whole burst and its count rate in 
the 8-20 keV band reaches about 10-40% of the maximum 
count rate in the same band during the GRB. 

As the study of the GINGA data shows evidence for a 
thermal origin, one can wonder whether these X-ray precur- 
sors are associated to the photospheric emission. This could 
be possible if the two following conditions are satisfied : 
-Condition (1) The ratio of the photospheric over internal 
shock count rate as defined by eq. ^ must be small in the 
gamma-ray range but greater than a few 10% in the X-ray 
range. The region of the fcTph /(l + z)-\/ / 7 plane where 
such a condition can be achieved (using the energy bands of 
the WATCH experiment) is shown in fig. We find that (i) 
the photospheric temperature kT p ^ + z) must lie in the 
X-ray band, which is easily obtained if A ~ 0.1. (ii) the ra- 
tio A/ / 7 must be above a minimal value which is typically 
about 0.1 and decreases when the peak energy E p increases. 
-Condition (2) As no activity (thermal or non thermal) is 
detected in the gamma-ray band during the X-ray precur- 



sor, the internal shock emission must start at the end of the 
precursor. We have shown in sec. 3.2 that the arrival time 
of the photons emitted by the layer ejected by the source 
at tinj when it becomes transparent can be approximated by 
t a ~ ti n j . The arrival time of photons emitted by the internal 
shocks due to the collisions between two layers emitted at 

(with r(i inj ) > r(4j)) is 



^inj ^inj ^ ^inj 



fa 



finj + ' 



t- ■ - t- 



(r(i inj )/r(*; nj )) -1 



(46) 



The only possibility to increase the delay between the 
beginning of the photospheric emission and the beginning 
of the internal shock emission is then to impose that the 
variability of the initial distribution of the Lorentz factor in 
the relativitic wind is initially low (F(fi n j)/r(t[ n j) — > 1) and 
increases during the wind production by the source. 

Condition (1) is easily achieved if the initial fraction of 
the energy released by the source under internal energy 
form is low. For instance, in fig. ^, one clearly sees a X-ray 
precursor lasting for about 10% of the total duration with 
an intensity of about 8% of the intensity at maximum (in 
the X-ray band). On the other hand, some of the precursors 
observed by GINGA and WATCH/GRAN AT have longer 
durations. This is where condition (2), which is probably a 
stronger constraint, is important. 

We propose the following interpretation for the pres- 
ence or absence of a precursor : it is necessary to have 
A <C 1 in order to suppress a too strong thermal gamma-ray 
emission from the photosphere. This naturally leads to a 
prompt thermal X-ray activity, which then could be very 
frequent in GRBs. However this activity is too weak to 
be easily detected when it occurs simultaneously to the 
bright non-thermal emission from the internal shocks. It 
is only when it appears as a precursor activity that it 
can be clearly identified. This can happen if by chance 
the relativistic wind is initially produced with a smooth 
distribution so that the internal shock activity is delayed. 
The expected features of such precursors are very close to 
the properties of the X-ray precursors observed by GINGA 
and WATCH/GRAN AT. 

To check the validity of this interpretation one clearly 
needs more precise detections of the X-ray prompt emis- 
sion of GRBs and especially a better characterization of the 
spectral properties of the X-ray precursors. If a black-body 
spectrum can be identified without any ambiguity, the cor- 
responding temperature will be measured, which would con- 
strain the A parameter. 

5.2 The optical photospheric emission 

To recover a dominant non-thermal gamma-ray emission 
we need A to be of a few percents or less. The correspond- 
ing Planck spectrum then peaks in th e X-ray band : for 
instance, the burst considered in sec 
temperature in the range kT p ^ ~ 
and fcTph ~ 0.6-30 keV for A = 0.01. It is interesting to 
estimate what is the photospheric emission for even lower 
values of A and if it could be dominant in the optical and 
produce a prompt optical flash comparable to that observed 



4.3 has a photospheric 
2-100 keV for A = 1 
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0.1 1 10 100 0.1 1 10 100 
kT ph /(l+z) [keV] kT ph /(l+z) [keV] 

Figure 7. In the fcT p h /(l + 2)— A// 7 plane, the stripped area 
shows the region where the photospheric thermal emission is neg- 
ligible in the 20-60 keV gamma-ray band of WATCH/GRANAT 
(count rate ratio smaller than 1% ) but can be detected in the 
8-20 keV X-ray band of the same experiment (count rate ratio 
greater than 10% ). Each panel corresponds to a particular value 
of the peak energy E p of the internal shock emission (respectively 
50, 200, 500 and 1000 keV). The solid, dotted and dashed curves 
respectively correspond to 2 = 0, 2 = 1 and 2 = 3. The other pa- 
rameters are Ljg = / 7 E = 10 51 erg/s, a = —1. and (3 = —2.25. 



in GRB 990123 (Akerlof et al. 1999) 



It is very unlikely that the photospheric emission peaks 
in the V band, because the photospheric temperature scales 
as A 1,/4 . To decrease fcT p h from 100 keV to 1 eV, a very 
unrealistic value of A = (10~ 5 ) 4 = 10~ 20 is required ! The V 
band then always lies in the Rayleigh- Jeans part of the pho- 

dn ph 



tospheric spectrum where ■ 



oc E. GRB 990123 has an av- 



eraged spectrum which is well reproduced by the Band func- 
tion with the following parameters: a — —0.6, j3 = —3.11, 
E p = 720 keV and photon fl ux 1.93 10~ 3 ph/s/cm 2 /keV 
at 1 MeV flBriggs et al. 1999| ). The redshift of the source 
is z — 1.6. The corresponding internal shock luminosity is 
Lis = /t E ~ 8.6 10 52 erg/s. If we now assume that the 
internal shocks have no other contribution in the optical 
range than that given by the Band spectrum, their flux in 
the V band (0.55 rim) is Fy S ~ 4.9 10~ 2 mjy, which is 
much too low to explain the optical flash reaching magni- 
tude m v ~ 9 (i.e. F v ~ 0.92 Jy ) observed by ROTSE. 
In the Rayleigh- Jeans regime, the corresponding flux due to 
the thermal photospheric emission is even lower : 



E 



pi, 



2.:! in ' I A 



3/4 



1/2 



/ fcT ph \ - 
U kevj 



mjy 



(47) 



We then find that the photospheric optical emission is much 
too weak to explain the ROTSE observations. This result 
is mainly due to the fact that the photospheric luminosity 
decreases much faster with A than the temperature. 



6 CONCLUSIONS 

In the framework of the internal shock model for gamma-ray 
bursts, we have computed in a detailed way the photo- 
spheric emission of an ultra-relativistic wind with a variable 
initial distribution of the Lorentz factor. We have compared 
the obtained spectrum and time profile to the non-thermal 
contribution of the internal shocks. Our main results are 
the following : 

(1) The photosphere in the standard fireball model is too 
hot and luminous. In the standard fireball model where the 
initial temperature of the fireball is about 1 MeV, the inter- 
nal energy is still large when the wind becomes transparent 
and the photosphere is therefore hot and luminous. The 
consequence is that the photospheric thermal component 
in the X-ray /gamma-ray range is in most cases at least as 
bright as the non-thermal component due to the internal 
shocks (even if the internal shock efficiency is high). This 
is in contradiction with the observations of BATSE and 
Beppo-SAX showing non-thermal spectra. 

(2) MHD winds are favored. Results in much better agree- 
ment with the observations are obtained when it is assumed 
that only a small fraction A of the energy released by the 
source is initially injected under internal energy form in a 
fireball. Most of the energy could for instance be initially 
under magnetic form, a large fraction of the Poynting flux 
being eventually converted into kinetic energy at large 
distances. For a typical internal shock efficiency of a few 
percents, values of A < 0.01 are required, which means that 
not more than 1% of the energy is initially deposited in the 
ejected matter (whose initial temperature is then of about 
a few hundreds keV). 

(3) X-ray thermal precursors can be obtained. A conse- 
quence of this strong assumption is that moderately low 
A (A ~ a few percents) lead to the presence of thermal 
X-ray precursors if the distribution of the Lorentz factor 
is not too variable in the initial phase of wind production. 
The characteristics of these precursors (spectral range, 
duration, intensity) are very comparable to the X-ray 
precursor activity observed in several GRBs by GINGA 
and WATCH/GRANAT. 

(4) The optical photospheric emission is very weak. For very 
small A values, the photospheric emission can be shifted 
to even lower energies. However, we have shown that it 
also becomes much too weak to explain the prompt optical 
emission observed by ROTSE in GRB 990123. 

A good test of the results presented in this paper would 
be the detection of X-ray precursors by an instrument with 
good spectral capabilities, so that a thermal origin could 
be firmly established. A determination of the photospheric 
temperature would put an interesting constraint on the A// 7 
ratio and then on the wind acceleration mechanism. More- 
over, if the photospheric thermal emission could be clearly 
detected (for instance in the soft X-ray range), it would pro- 
vide a direct information about the initial distribution of the 
Lorentz factor in the wind before the internal shocks start. 



ACKNOLEDGMENTS 

F.D. acknowledges financial support from a postdoctoral fel- 
lowship from the French Spatial Agency (CNES). 



The expected thermal precursors of GRBs 11 



REFERENCES 

Abramowicz M. A., Novikov I. D., Paczynski B., 1991, ApJ, 
369, 175 

Akerlof C, Balsano R., Barthelemy S., Bloch J., Butter- 
worth P., Casperson D., Cline T., Fletcher S., Frontera 

F. , Gislcr G., Heise J., Hills J., Kehoe R., Lee B., Mar- 
shall S., McKay T., Miller R., Piro L., Priedhorsky W., 
Szymanski J., Wren J., 1999, Nature, 398, 400 

Band D., Matteson J., Ford L., Schaefer B., Palmer D., 
Teegarden B., Cline T., Briggs M., Paciesas W., Pendleton 

G. , Fishman G., Kouveliotou C, Meegan C, Wilson R., 
Lestrade P., 1993, ApJ, 413, 281 

Baring M. G., Harding A. K., 1997, ApJ, 491, 663 
Blandford R. D., McKee C. F., 1976, Physics of Fluids, 19, 
1130 

Blandford R. D., Znajek R. L., 1977, MNRAS, 179, 433 
Briggs M. S., Band D. L., Kippen R. M., Preece R. D., 
Kouveliotou C, van Paradijs J., Share G. H., Murphy 
R. J., Matz S. M., Connors A., Winkler C, McConnell 
M. L., Ryan J. M., Williams O. R., Young C. A., Dingus 
B., Catelli J. R., Wijers R. A. M. J., 1999, ApJ, 524, 82 
Chevalier R. A., Li Z., 2000, ApJ, 536, 195 
Daigne F., Mochkovitch R., 1998, MNRAS, 296, 275 
Daigne F., Mochkovitch R., 2000, A&A, 358, 1157 
Djorgovski S. G., Kulkarni S. R., Bloom J. S., Frail D. A., 
Harrison F. A., Galama T. J., Reichart D., Castro S. M., 
Fox D., Sari R., Berger E., Price P., Yost S., Goodrich R., 
Chaffee F., 2001, in Gamma-ray Bursts in the Afterglow 
Era The GRB Host Galaxies and Redshifts. p. 218 
Frail D. A., Kulkarni S. R., Nicastro S. R., Feroci M., Tay- 
lor G. B., 1997, Nature, 389, 261 
Frail D. A., Kulkarni S. R., Sari R., Djorgovski S. G., Bloom 
J. S., Galama T. J., Reichart D. E., Berger E., Harrison 
F. A., Price P. A., Yost S. A., Diercks A., Goodrich R. W., 
Chaffee F., 2001, ApJ, 562, L55 
Frontera F., Amati L., Costa E., Muller J. M., Pian E., Piro 
L., Soffitta P., Tavani M., Castro- Tirado A., Dal Fiume 

D. , Feroci M., Heise J., Masetti N., Nicastro L., Orlandini 
M., Palazzi E., Sari R., 2000, ApJS, 127, 59 

Frontera F., Amati L., Vietri M., Zand J. J. M. L, Costa 

E. , Feroci M., Heise J., Masetti N., Nicastro L., Orlandini 
M., Palazzi E., Pian E., Piro L., Soffitta P., 2001, ApJ, 
550, L47 

Ghisellini G., Celotti A., Lazzati D., 2000, MNRAS, 313, 
LI 

Kobayashi S., Piran T., Sari R., 1997, ApJ, 490, 92 
Lithwick Y., Sari R., 2001, ApJ, 555, 540 
Lyutikov M., Usov V. V., 2000, ApJ, 543, L129 
Meszaros P., Rees M. J., 2000, ApJ, 530, 292 
MacFadyen A. L, Woosley S. E., 1999, ApJ, 524, 262 
Meszaros P., Rees M. J., 1997, ApJ, 482, L29 
Mochkovitch R., Hernanz M., fsern J., Martin X., 1993, 

Nature, 361, 236 
Murakami T., Inoue H., Nishimura J., van Paradijs J., Fen- 

imore E. E., 1991, Nature, 350, 592 
Narayan R., Paczynski B., Piran T., 1992, ApJ, 395, L83 
Paczynski B., 1998, ApJ, 494, L45 
Panaitescu A., Kumar P., 2001, ApJ, 554, 667 
Piran T., 1999, Physics Reports, 314, 575 
Preece R. D., Briggs M. S., Mallozzi R. S., Pendleton G. N., 

Paciesas W. S., Band D. L., 2000, ApJS, 126, 19 



Ramirez-Ruiz E., MacFadyen A. L, Lazzati D., 2002, MN- 
RAS, 331, 197 

Rees M. J., Meszaros P., 1994, ApJ, 430, L93 

Rhoads J. E., 1999, ApJ, 525, 737 

Sari R., Piran T., Narayan R., 1998, ApJ, 497, L17 

Sazonov S. Y., Sunyaev R. A., Terckhov O. V., Lund N., 
Brandt S., Castro- Tirado A. J., 1998, A&ASS, 129, 1 

Spruit H. C, Daigne F., Drenkhahn G., 2001, A&A, 369, 
694 

Thompson C, 1994, MNRAS, 270, 480 

van Paradijs J., Groot P. J., Galama T., Kouveliotou C, 
Strom R. C, Telting J., Rutten R. G. M., Fishman G. J., 
Meegan C. A., Pettini M., Tanvir N., Bloom J., Pedersen 
H., Nordgaard-Nielsen H. U., Linden- Vornle M., Melnick 
J., van der Steene G., Bremer M., Naber R., Heise J., 
in 't Zand J., Costa E., Feroci M., Piro L., Frontera F., 
Zavattini G., Nicastro L., Palazzi E., Bennet K., Hanlon 
L., Parmar A., 1997, Nature, 386, 686 

Wijers R. A. M. J., Rees M. J., Meszaros P., 1997, MNRAS, 
288, L51 

Woosley S. E., 1993, ApJ, 405, 273 



